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. ABSTRACT. We show that a conditional characteristic function of gener- 

r-j ! alized Levy stochastic areas can be viewed as a product a selfdecomposable 

distribution (i.e., Levy class L distribution) and its background driving char- 
acteristic function. This provides a stochastic interpretation for a ratio of 
some Bessel functions as well as examples of characteristic functions from 
van Dantzig class. 
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ov 

O ' 1- An introduction. For a planar Brownian motion B t = (Z t ,Z t ) and a 

stochastic area process 



A u = I Z s dZ s - Z s dZ s , w > 0, 

Paul Levy (1950) proved that its conditional characteristic function is of the 
form 



tu \(x\^ 

E\e ltAu \B u = a] = expf- — (tu cothtu - 1)1, t e R, 

L sinhtw FL 2u y n 

where a£l 2 and u > are fixed. Thus, in particular, 

E[e ltAu \B u = Vu)) = — : exp\-(tu cothtu - 1)1, teR. (1) 

smh tu 



*This work was completed while the author was visiting (2002/03) Wayne State Uni- 
versity, Detroit, USA. 
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For another proof of Levy stochastic area formula and related topics involving 
more stochastic argument we refer to Williams (1976) and Yor (1980). 

In Jurek (2001), p. 248, it was showed that in the formula (1) one has a 
product of the class L characteristic function btj sinh bt, {b G K is a fixed pa- 
rameter), and a characteristic function exp[— (bt coth bt — 1)] of its background 
driving random variable that is from the class hi of s-selfdecomposable dis- 
tributions. Consequently, by Proposition 3, in Iksanov, Jurek and Schreiber 
(2003), we have that (1) is again in class L. In this note we prove that a 
similar result holds for a generalized Levy stochastic areas; cf. Theorem 1. 
Here hyperbolic trigonometric functions are replaced by the modified Bessel 
function of first kind l v . As a by-product we get a formula for a ratio of such 
functions and provide examples of van Dantzig class V; cf. Proposition 1, 
Corollaries 2-4 and the formula (12). 

Bessel functions I u (t), as completely monotone functions in \ft or ^/v or 
as functions of a pair (t, u), appeared in many areas of probability and can be 
interpreted by some stochastic functionals; cf. Ismail and Kelker (1979), but 
in particular, cf. Pitman and Yor (1981) and references there- in. In this note 
we provide stochastic way of looking at I u {t), and their ratios, as function in 



Our proofs here rely on different representations of Bessel functions and 
therefore for ease of reading, in the Appendix, we have collected all needed 
formulas. This also makes this paper more self-contained. 

2. Selfdecomposable and s-selfdecomposable distributions. 

Self decomposable random variables X (or their probability distribution 
or their characteristic functions 0) are defined as limits of normalized partial 
sums of infinitesimal independent random variables. Thus they are infinitely 
divisible, in short: X G ID. Let L stand for the class of all either selfde- 
composable random variables X or their probability distributions ji or their 
characteristic functions 0. Then 



and this is so called the selfdecomposability property. One of the major char- 
acterizations of the class L, called the random integral representation is as 
follows 



t G R. 



<p G L iff V(0 < c < l)3(char.f. </>(■) = <j>(c-)ip c (-) 



(2) 




(0,oo) 



e- s dY(s), with £[log(l + |y(l)|)] < oo, 



(3) 
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where = means equality in distribution and Y(-) is a Levy process called 
background driving Levy process; in short: Y is BDLP of X or X is driven 
by Y; cf. Jurek and Mason (1993), Chapter 3. For BDLP identification 
purposes note that if ip is the characteristic function of Y{1) in (3) then 



In particular, class L characteristic functions are differentiable for t ^ and 
if) is infinitely divisible characteristic function with a finite logarithmic mo- 
ment; cf. Jurek (2001), Proposition 3. Finally let us recall that (L, *) is a 
closed topological subsemigroup of the semigroup ID. 

Class L is quite rich and contains among others: all stable laws, gamma 
j a ,x, Laplace (double exponential), t-Student, F-Fisher, log-normal, log-F, 
log 1 1 1 , inverse Gaussian, Barndorff-Nielsen generalized hyperbolic, general- 
ized gamma, etc. 

Another class of infinitely divisible distributions needed here is the class 
U of s-selfdecomposable probability measures. One of the descriptions of U is 
as follows 

ijj G U iff V(0 < c < l)3(char.f. p c G ID) ip(t) = ^ c (ct)p c (t), t G R, (5) 

where the power is well defined as ip(t) ^ 0. To the above we refer to as 
the s-selfdecomposability probability. From (5), in particular, we infer that 
(U, *) is also a closed convolution topological semigroup. In fact, we have 
the inclusions IcWc ID. Similarly, for s-selfdecomposable distributions 
we have analogous random integral representation 



and as before we refer to Y as the BDLP for s-selfdecomposable X. From 
(6) we infer that if and ip are characteristic functions of X and Y(l), 
respectively, then 

if>(t) = exp[(flog0(t))'], t^O and ^(0) = 1. 

For those and other relations between the classes L and U we refer to 
Jurek (1985) or Iksanov, Jurek and Schreiber (2003). 



ip(t) = exp[t(log0(t))'], * ^ 0, and ^(0) = 1. 



(4) 




(6) 
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3. Stochastic areas and selfdecomposability. For p > —1/2, let 
us define a Gaussian process V p and the corresponding generalized Levy 
stochastic area as follows 



V? = (Vf, Vf) := t-' p f dB s , and A p := /' V? dV? - V* dV* 

Jo Jo 

where B s , s > 0, is the planar Brownian motion. In Biane and Yor (1987) 
(see also Yor (1989), p. 3052 or Duplantier (1989)) the following generaliza- 
tion of Levy stochastic area formula was proved 

1 |Vl - aJ -2^+l)/,(|A|) eXPi 2 |A| mx\) 1 {7) 

where I v (z) is the modified Bessel function of the index v := p+1/2 > and 
a G M 2 is fixed. In the 'products' in (7) we recognize a similar formula as in 
(1) and , moreover, the components have analogous interpretations. Namely 
we have the following result. 

THEOREM 1. (%) For v > -I, functions 

f 

BJt) := ; — , tsl, are class L characteristic functions. 

w 2 v F(u + 1) i„(ty 

More explicitly, these are the characteristic functions of random variables of 
the form 

oo 

X u := z~l rjk, where i] k are i. i. d. Laplace rv's, 
k=i 

and z„ t k, k = 1,2, are the zeros of the function z~ u J v (z). Here J u is the 
Bessel function of the first kind. Furthermore 

p oo 

logB v (t) = - (1 - costx)(J2 e-^ xl )/\x\dx. 
Jr\{o} k=1 

(ii) Functions 

K(t) := exp [ — t -r - ] , t G R, are class U characteristic functions. 
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Moreover, these are the background driving characteristic functions for B v , 
i.e., log&j,(i) = td(\ogB u (t))/dt, t ^ 0. Furthermore, 



it/2) I u+1 (t) = I u (t) / (1 - coste)( V \z v , k \e-\*"*\ x )dx. 
Jo k=i 



(8) 



(Hi) For all positive constants a and (3, the functions B«(t) ■ &£(*) are 
s-self decomposable (class U distributions). Moreover, for < /3 < a the 
resulting characteristic functions are self decomposable (class L characteristic 
functions). 

(iv) For v > —1/2, functions \/B v (it) = r ^ 1 it^ 1 i G 1, are charac- 
teristic functions corresponding to the probability densities 

These are not infinitely divisible characteristic functions. 
Proof. From A. 3 (a) in the Appendix, we have 

where z v ^ are real zeros of the function z~ u J u (z), ordered accordingly to 
their absolute values; cf. also A. 1(a) in the Appendix. Since Laplace (also 
called the double exponential) random variable r\ has characteristic function 



/oo r -\ x \ 

e *tx 2 -i e -\x\ dx = ( 1 + t 2yi = exp / ( cog fx _ l) dx, 
■00 Jr\w\ fI 



cos trr — 1) — 

lR\{0} ' 

we easily conclude that it is class L distribution. Furthermore, 

2 - 00 



= II + — ) _i = CX P / (costs- l)(J^ e 

fc=l M{o} fc= i 



-^■*H x l)/|a;| (10) 



cf. Jurek (1996), p. 177 for more details. Since class L is closed under dila- 
tions (multiplication of rv by scalars) and the weak convergence therefore X u 
has class L distribution as well, with the characteristic function B v (t), which 
proves the part (i). 
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Part (ii). By Proposition 3 in Jurek (2001) or see (4), for any class L 
characteristic function <f>(t) we have that ip(t) := exp(td(log<f)(t)) / dt) exists 
(for t ^ 0) and it is the characteristic function of BDLP. Since z ^ Iu(z) — 
vl v (z) = zl v+ i(z), by A. 8, we get 

tdQogB v (t))/dt = u-t!j& = -t = log 6,(0, 

which proves that b v (t) are indeed characteristic functions of the BDLP for 

rv X v . Using (4) and Jurek (1996), we get the formula (8). Finally, since X v 

are given by series on independent Laplace rv therefore by Proposition 3 in 

Iksanov, Jurek and Schreiber (2003) to infer that b v (t) is from the class U of 

s-selfdecomposable distributions. 

Part (iii). Since L dU and both classes are closed under taking positive 

powers therefore B$(t) ■ &£(t) G U, which proves the first part of (iii). Let 

a > f3 > 0. Since b v {t) G U therefore B u (t) ■ b u (t) G L, by Theorem 1 in 

— — i 

Iksanov, Jurek and Schreiber (2003). Since also By (t) G L and L is a 
semigroup, we conclude that B v (t) a ^ ■ b v (t) G L. Taking powers we get the 
second part of (iii). 

For part (iv), using A. 5 and part (i) of the theorem, we get 

—L- = — r -^ ^ f 1 e itx (l - a?)"-vdx, (11) 

B„(it) B{\,u+\)J^ K ' K } 

where B(x,y) = is Euler's beta function (Euler's integral of the first 

kind). Since, by A.l, J u (t) has real zeros therefore it can not be infinite 
divisible. This completes the proof of the Theorem. 

For some related information on random variables of the form X v and the 
corresponding Dirichlet series we refer to Jurek (2000). Here, in particular, 
we have the following conditions on zeros of Bessel functions. 

COROLLARY 1. For zeros z v ^, k = 1,2, of the function z~ v J v (z) with 
v > —1, we have 

y~] \z vk \ / log \x\ e~\ z "> k ^dx < oo, 

k ' J{\x\>l) 

J2\z uk \ [ log(l + \x\ 2 ) e-\ z ^dx < oo. 

k ' JQx\>0) 
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[Both integrals can not be expressed in simple functions. The first involves 
the incomplete gamma function and the second sine and cosine integrals; cf. 
Gradshetyn and Ryzhik (1965) formulae 4.358 and 4.338, respectively.] 

Proof. From (i) in Theorem 1 we infer that X v has Levy spectral M of 
the form 

oo „ 

M(A) = V / e -1 *"'* 11 * 1 \x\~ l dx, for all Borel subsets of K \ {0}. 
k=i Ja 

And from random integral representation (3) or from (8) in Theorem l(ii) 
we have 

M(A) = / N(e s A)ds and N(A) \zu,k\e~ MM dx, 

</(0,oo) k=l^ A 



and A is a Levy spectral measure of BDLP Y with finite logarithmic moment 
on {\x\ > 1}. This gives the first condition in Corollary. For the second one 
we use the fact that 

^— 2 M{dx) < oo, 
{M>o} 1 + x 

and the integration by parts formula. 

REMARK 1. Biane and Yor (1987) expanded Brownian path along Legendre 
polynomials and computed conditional characteristic function of stochastic 
area given values of few first coefficients. In that case again we have formulas 
involving the products (7); cf. formula (4.10) in Yor (1989). 

REMARK 2. The fractions I u+ i/I u of Bessel functions appear in the formula 
for the background driving characteristic functions b v for the selfdecompos- 
able characteristic functions B v (t); cf. Theorem l(ii). In Jurek (2001), in 
Example 2, it was found that functions 

t-^exp[-|t|^M], teR, 

involving fractions of modified Bessel functions K u , are background driving 
characteristic functions as well. This time they correspond to Student's t- 
distributions with 2u degrees of freedom. 
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4. van Dantzig distributions and selfdecomposability. 

Recall that an analytic characteristic function <p(z) defined in a strip —a < 
Qz < b, for some a, b, > 0, belongs to van Danztig class V, if ip(t) := l/<p(it) 
can be extended to a characteristic function; cf. Lukacs (1968). Equivalently, 

i/j(t) ■ <p{it) = 1, tel, (van Dantzig pair (0,^)). (12) 

If ip(t) = <j){t) we say is self-reciprocal characteristic function in van Danzitg 
class V. The three elementary examples of pairs (ip, 0) are: 

(cost, (cosht)- 1 ), and (e" rf2 ,e" c * 2 ) with c > 0. 

The normal characteristic function is an example of self-reciprocal one. Here 
is another example. 

PROPOSITION 1. For v > -1/2, the ratios t € R, are non- 

infinitely- divisible self -reciprocal characteristic functions from van Danztig 
class V. 

Proof. From Theorem 1, parts (i) and (iv), we have that the ratio is indeed 
a characteristic function. Since it has real zeros therefore ^4rv can not be 
infinitely divisible. Its self-reciprocal property follows from the identities A. 2 
in the Appendix and the fact that B u are even functions; cf. in the Appendix 
formulae A. 3 or A. 3a. Thus the proof is complete. 

COROLLARY 2. All characteristic functions B v , for v > —1/2, belong to 
van Dantzig class V. 

This is a consequence of Theorem 1, parts (i) and (iv). Here are more 
explicit formulae of class L and van Dantzig characteristic functions. 

COROLLARY 3. For n = 0, 1, 2, ... and t e R, functions 

((2n + l)!!)- 1 t 2n+1 



B n+h {t) 



{E{k:n-keven} a k,nt n k ) SUlll t + {E{k: n-k odd} a k,nt n k ) COsh t 



where a kj , a := ( — 2 kk\{n-k)\ ' ® — ^ — n > are (characteristic functions of class 
L distributions that belong to van Dantzig class V. 
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The proof is a combination of the formulae 




from A. 6a and T(n + |) = < - 2 " 2ra 1 ^ !! , n > 1, and = y^, with the part 
(i) of Theorem 1. 

COROLLARY 4. Forn = 0, 1, 2, ... and van Dantzig class V characteristic 
functions B n+ i{t), we have the equality 



This is a consequence of the identity I n+ i(z) = z n+ ? \j\ i^zY (~T £ )' c ^ 



A. 7a, in the Appendix, and the Theorem l(i). 

REMARK 3. Note that our class L characteristic functions Bi(t), Bs(t) and 
Bh (t) are exactly the examples (of the van Dantzig class T>) given in Lukacs 
(1968), on p. 121, formulae (14. 2a), (14. 3a) and (14.4a). However, they were 
obtained there by different inductive procedures. 

REMARK 4. It might be worthy to notice that Lukacs' 1968 examples of 
class V characteristic functions, given below his formula (13) on p. 120, 
are class L characteristic functions. The same applies to those given at the 
bottom of page 121 (convolutions of normal distributions with several gamma 
distributions). All in all it might be true that in each non self- reciprocal pair 
((f), ip) of characteristic functions satisfying van Dantzig condition (12) one 
of them is selfdecomposable. Or if it is not the case, then characterize that 
part of V that possesses this property. 

REMARK 5. In search for examples of van Danztig pairs one may follow 
recent paper by de Meyer, Roynette, Vallois and Yor (2002). In Theorem 4.1 
Authors have the equality 



where 



(tdt) 



d' 



,71 
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for a standard Brownian motion B t ,£ > 0, and some independent of it stop- 
ping times T. 

Acknowledgement. Author would like to thank Professor Marc Yor, 
from University of Paris VI, for some valuable comments and for providing 
the references [14], [16] and [17]. 

5. Appendix. Recall that the Bessel functions are the solutions of the 
differential equation 

^ + i^ + (l--)^ = 0. 

dz 2 z dz z 2 v 

Our main interest here is in the Bessel function of the first kind J u (z) and 
their modified versions I u (z), with z e C. 

For ease of reading we collect here some formulas for Bessel functions. 
They are either taken directly from Gradshteyn and Ryzhik (1965), and in 
that case they are labelled by bold face numbers, or are combination of those. 

A.l^) = gn:,(l-0 ^-1,-2,-3,..., 8.544 

where z v>m are zeros of the function z~ v J u (z), ordered 
by absolute value of their real parts. 

A. la For v > — 1 all zeros z v<m are real numbers. 8.541 

A. 2 I v (z) = e~^ ih 'J^e^z), —n < argz < n/2, and 

I v \z) = e^J^e'^z), tt/2 < arg^ < 7T. 8.406 

A.3 /^) = ^n~= 1 (l + £), ^-1,-2,-3,..., 

This is derived from the factorization A.l with a usage of A. 2. 

A.3a I v (z) = (z/2Y Er=o m^kWy 8 ' 445 

AA J ^ z) = roSfe) /i e * 2<(1 - t2)v ^ dt > > "5> 8 - 41 ( 10 ) 

A.5 Uz) = Y ^ W) ! l _ l e^\l~t 2 r-Ut, &/>-§, 8.431(1) 

A - 6 WfcW^ElU*^^ 8 ' 467 
with dk, n '■= (~l) fc 2 fc fc!(n-fc)! ' < A; < n, which may be rewritten as 
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A-6a I ±{n+l) {z) = y/lz-^) 

■ [{J2{k:n-keven} a k,nZ n - k ) Smh Z + ( J2 {k :n-kodd} a k,n Z n ~ k ) COsh z] 

\ ^d z y ^ s the n -times composition of an operator / — > | /'(z)]- 

A.7a/ B+ i(z) = ^ > /|^(^) 

It is a consequence of the formula A. 7 and A. 2. 

A.8 z£l v (z)-ul v (z)=zl v+1 (z) 8.486(4) 
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